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Since Dicke’s seminal paper on coherence in spontaneous radiation by atomic ensembles, superradiance has 
been extensively studied. Subradiance, on the contrary, has remained elusive, mainly because subradiant states 
are weakly coupled to the environment and are very sensitive to nonradiative decoherence processes. Here we 
report the experimental observation of subradiance in an extended and dilute cold-atom sample containing a 
large number of particles. We use a far detuned laser to avoid multiple scattering and observe the temporal 
decay after a sudden switch-off of the laser beam. After the fast decay of most of the fluorescence, we detect a 
very slow decay, with time constants as long as 100 times the natural lifetime of the excited state of individual 
atoms. This subradiant time constant scales linearly with the cooperativity parameter, corresponding to 
the on-resonance optical depth of the sample, and is independent of the laser detuning, as expected from a 
coupled-dipole model. 


Despite its many applications, ranging from astro- 
physics [U to mesoscopic physics mia and quantum in¬ 
formation technology [1], light interacting with a large 
ensemble of N scatterers still bears many surprising 
features and is at the focus of intense research. For 
N = 2 atoms placed close together, the in-phase oscilla¬ 
tion of the induced dipoles produces a large, superradiant 
dipole, whereas the out-of-phase oscillation corresponds 
to a subradiant quadrupole. Generalizing for V ^ 2, 
Dicke has shown that, for samples of a small size com¬ 
pared to the wavelength of the atomic transition, the 
symmetric superposition of atomic states induces super¬ 
radiant emission, scaling with the number of particles N, 
whereas the antisymmetric superpositions are decoupled 
from the environment, with vanishing emission rates, 
which corresponds to subradiance [5]. 

Dicke superradiance has been extensively studied in 
the 1970s [BHS] but the observation of its counterpart, 
subradiance, has been restricted to indirect evidence of 
modified decay rates in one particular direction or in 
systems of two particles at very short distance [TOHT^ . 
One challenge for the observation of subradiance by a 
large number of particles is the fragile nature of these 
states, which require protection from any local nonra¬ 
diative decay mechanism m- Furthermore, contrary to 
the two-atom case, for which the distance between atoms 
has to be small compared to the wavelength, for V 3> 2, 
the retarded, long-range resonant dipole-dipole interac¬ 
tion [14] gives rise to super- and subradiant effects (“co¬ 
operative scattering”) also in dilute samples, with inter¬ 
atomic distances much larger than the wavelength, and 
corresponding large system sizes. Since, for N > 2, the 
Hamiltonians for short and long-range interactions do 
not commute, the collective eigenstates due to the long- 
range interactions are suppressed by short-range interac¬ 
tions [Sj. These short-range or near-field effects (or ‘van 
der Waals dephasing’) thus need to be avoided in this 
case. As a consequence - and maybe counterintuitively 


- a large and dilute sample of interacting dipoles is the 
most appropriate system for the observation of A-body 
subradiance. 

In this regime, and in the weak excitation limit 
(“single-photon superradiance” ) dSHni, it has been 
shown that the superradiant enhancement of the emis¬ 
sion rate scales as the cooperativity parameter N/AI, 
where M is the number of available modes for the electro¬ 
magnetic radiation [IBHU- For a spherical sample of ra¬ 
dius R, M ^ {koR)^, where ko = 27r/A, this cooperativ¬ 
ity parameter is proportional to the peak on-resonant op¬ 
tical depth of the atomic cloud, given by bo = 3N/{koR)‘^ 
for a cold-atom cloud with a Gaussian density distribu¬ 
tion of rms radius R. This number can be large even 
at low density. In a recent work [52], we used a coupled- 
dipole model to generalize this result to subradiance (see 
also the Supplemental Material [23]). In this Letter, we 
report the experimental observation of subradiance in 
this weak-excitation, dilute- and extended-sample limit. 

In our experiment, we load A Ri 10® ®^Rb atoms from 
a background vapor into a magneto-optical trap (MOT) 
for 50 ms. A compressed MOT (30 ms) period allows for 
an increased and smooth spatial density with a Gaus¬ 
sian distribution of rms size i? Ri 1 mm (typical density 
p pa 10^^ cm“®) and a reduced temperature T « 50 p,K. 
We then switch off the MOT trapping beams and mag¬ 
netic field gradient and allow for 3 ms of free expansion, 
used to optically pump all atoms into the upper hyper- 
fine ground state F = 2. Next, we apply a series of 12 
pulses of a weak probe beam (waist w = 5.7 mm), lin¬ 
early polarized and detuned by d = (w — a;o)/r from the 
closed atomic transition F = 2 —>■ F' = 3. Here uj is the 
frequency of the laser, ojq the frequency of the atomic 
transition (of wavelength A = 27rc/a;o = 780.24 nm) 
and r/27r = 6.07 MHz its linewidth. Note that when 
we varied the detuning, we also varied the laser inten¬ 
sity accordingly in order to keep the saturation param¬ 
eter constant at s ~ 4.5 x 10“^. The pulses of duration 
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Fig. 1. (color online). Principle of the experiment. A 
large probe laser illuminates the atomic sample for 30 ^s 
and is switched off rapidly. The fluorescence at ~ 35° 
is collected by a hybrid photomultiplier (HPM) and 
recorded on a multichannel scaler (MCS). The experi¬ 
ment is repeated 500000 times. At each cycle, 12 pulses 
are recorded during the free expansion of the cloud, al¬ 
lowing the on-resonance optical depth to vary. 


30 fis and separated by 1 ms are obtained by using two 
acousto-optical modulators in series to reach an extinc¬ 
tion ratio better than 10“"^. The 90%-10% fall time at 
the switch-off is ^ 15 ns, limiting the possibility of study¬ 
ing superradiance, but convenient for detecting subradi¬ 
ance. Between subsequent pulses of each series, the size 
of the cloud increases because of thermal expansion, and 
the atom number decreases due to off-resonant optical 
pumping into the F = I hyperfine state during each 
pulse. The corresponding change of the on-resonant op¬ 
tical depth bo allows us to conveniently measure the de¬ 
cay of the fluorescence as a function of bo and investigate 
whether bo is the relevant scaling parameter [23]. After 
this stage of expansion and measurement, the MOT is 
switched on again and most of the atoms are recaptured. 
The complete cycle is thus short enough to allow the sig¬ 
nal integration over a large number of cycles, typically 
^ 500000 (complete acquisition time ^ 14h per run). 
The scattering of the probe beam is collected by a lens 
with a solid angle of 5 x 10“^ sr at 0 « 35° from the 
incident direction of the laser beam (see Fig. 1). We use a 
hybrid photomultiplier (Hamamatsu HPM R10467U-50) 
in the photon-counting regime, without any measurable 
amount of afterpulsing, which would considerably mask 
signatures of subradiance. The signal is then recorded on 
a multichannel scaler (MCS6A by FAST ComTec) with 
a time bin of 1.6 ns, averaging over the cycles. The co- 
operativity parameter bo corresponding to each pulse is 
calibrated by an independent measurement of the atom 
number, cloud size and temperature using absorption 
imaging [23]. 

Typical data are shown in Fig. 2. The signal is normal¬ 
ized to the steady-state fluorescence level and we focus 
on the switch-off period to highlight the slow fluores¬ 
cence decay. In Fig. 2(a), the detuning S of the probe 
beam is kept constant and the different decay curves 


correspond to different values of bo, obtained in a sin¬ 
gle run. On the contrary, Fig. 2(b) shows data taken 
with different detunings but for the same bo. In both 
cases, most of the fluorescence decays fast (note the log¬ 
arithmic scale of the vertical axis), but a slow decay is 
clearly seen well above the noise floor (slightly below 
10“^). We stress that fluorescence can be detected at 
very large delays, as can be seen from the time axis, in 
units of Tat = F”^ = 26 ns. We attribute this slow decay 
to subradiance in the single-photon (or weak excitation) 
regime, as predicted in Ref. [22]. 

A qualitative analysis of the two figures clearly shows 
different behaviors. As bo is varied, the slow decay rate 
changes, whereas its relative amplitude stays approxi¬ 
mately the same. The exact opposite happens when we 
change the detuning, keeping bo fixed. For a quantita¬ 
tive analysis, we fit the slow tail at long delays by an 
exponential decay with two free parameters: the time 
constant Tsub and its relative amplitude Agub [23] . We 
systematically studied how these parameters depend on 
bo and S. The result of this analysis is presented in Fig. 3. 
In Fig. 3(a) we plot the subradiant time constant as a 
function of bo for different detunings. The collapse of 
all points on the same curve clearly indicates that the 
slow decay rate does not depend on the detuning [see 
also Fig. 3(b)]. This demonstrates that this slow decay 
is not a multiple-scattering effect, such as previously ob¬ 
served radiation trapping |24j , which depends on the op¬ 
tical depth at the laser frequency b(S) oc bo/(l -I- 4(5^), 
with a strong dependance with <5. The second feature in 
Fig. 3(a) is the linear increase of Tgub with bo, up to time 
constants as long as Tgub lOOrat. This is perfectly con¬ 
sistent with the predictions of the coupled-dipole model 
for subradiance [23]. We note that for large negative de¬ 
tunings, one has to take into account the variation of bo 
during the pulse series induced by the cloud expansion 
together with a significant contribution of atom losses by 
off-resonant hyperfine pumping. This allowed us to test 
different combinations of IV and F as scaling parame¬ 
ters (see Fig. ?? in the Supplemental Material [23]). The 
comparison showed that the best collapse has been ob¬ 
tained with N/R^ (X. bo &s the scaling parameter, which 
demonstrates that bo is indeed the relevant cooperativity 
parameter in the regime of a dilute and extended sam¬ 
ple, as expected from the ratio of the number of atoms to 
the number of available electromagnetic modes radiating 
from the sample. Finally, we show in Fig. 3(c) the rela¬ 
tive amplitude Agub of the slow decay. As was already 
seen in Fig. 2(b), this amplitude is much larger near res¬ 
onance, and seems to reach a plateau for large detunings. 
We have checked that subradiant decay is still visible at 
a larger detuning, up to d = —11 [23]. This is in line 
with the coupled-dipole model, in which the weight of 
the long-lived modes are enhanced near resonance and 
the weight of all of the collective modes becomes inde¬ 
pendent of the detuning at large detuning. 

As a long lifetime can also occur due to multiple scat¬ 
tering, when the optical depth b{5) ^ 1 [53], we inves- 
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Fig. 2. (color online). Slow decay of the fluorescence after switching off the probe laser. The vertical scale is normalized 
by the steady-state detected power P(0) and the time scale is normalized by the atomic lifetime of the excited state 
Tat. Without any collective effect (single-atom physics), the decay would be given by P{t) = P(0) exp(—f/xat) (the 
black dotted line), (a) Several data are shown for different on-resonance optical depths bo and the same detuning 
S = —6 (in units of F). The time constant increases with bo- (b) Several data are shown for different detunings and 
the same bo = 108 ± 5. The time constant remains unchanged, but the relative amplitude of the subradiant decay 
decreases as the detuning increases. 


tigated the decay time close to the atomic resonance 
- to study how subradiance compares with radiation 
trapping. In the range of the accessible experimental 
values [Fig. 3(b)], no marked difference of the decay 
times around resonance is visible, even though a small 
difference was predicted in Ref. [52] . The interplay of ra¬ 
diation trapping and subradiance near resonance is still 
an open question and will be the subject of our further 
experiments. 

We have also studied the effect of the probe intensity 
and checked that, at low saturation parameter, the ob¬ 
served subradiance is independent of the intensity (see 
Fig. |6| in the Supplemental Material [23]b which vali¬ 
dates the use of the coupled-dipole model in the weak- 
excitation limit. We finally also excluded the possibility 
that residual near-resonant light might always be present 
and induce a slow decay due to radiation trapping, thus 
mimicking off-resonant subradiance |23j . 

To summarize, we have presented the first direct sig¬ 
natures of subradiance in a large system of resonant 
scatterers. We have shown that in the regime of dilute 
and extended samples, the subradiant decay rate is gov¬ 
erned by a cooperativity parameter defined as the ratio 
of the number of scatterers and the sample size squared, 
which conveniently corresponds to the on-resonance op¬ 
tical depth. This observation of subradiance opens in¬ 
teresting questions, including the robustness of subradi¬ 
ance against decoherence mechanisms or the possibility 
of controlling the population of the subradiant modes by 
an appropriate temporal or spatial shaping of the driv¬ 
ing laser or of the atomic levels. If the subradiant states 
can be manipulated with sufficient control [25] , their iso¬ 
lation from the environment might be exploited as a 
resource for quantum information or quantum metrol¬ 
ogy [26]. As subradiance goes hand in hand with super¬ 


radiance, simultaneous recording of fast and slow decays 
would be a beautiful illustration of the cooperative scat¬ 
tering envisioned by Dicke. By using a stronger laser 
drive, it would also be possible to access a larger sub¬ 
space of the full Hilbert space, addressing the possibility 
of a photon-blockade effect m- 

In addition to quantum optics, our observation is also 
relevant to mesoscopic physics [28l|29], a community 
less familiar with Dicke physics. One major challenge 
in this field is the observation of strong localization of 
light, the analogy for classical waves of Anderson local¬ 
ization of electrons 120]. Previous experimental obser¬ 
vations of a decay of scattered light slower than pre¬ 
dicted by the diffusion equation have been used as a 
signature of Anderson localization in dielectrics [311132] . 
Our results show that it cannot always be the case with 
cold atoms. Similarly, recent numerical simulations con¬ 
sidering point-dipole resonant scatterers study the col¬ 
lective modes of the effective Hamiltonian of the sys¬ 
tem and, in particular, their lifetimes [331137] . Our work 
shows that Dicke subradiance can also be at the origin of 
very long lifetimes and that a careful analysis is required 
to distinguish subradiant from localized modes [37]. Fi¬ 
nally, the combination of subradiance with disorder act¬ 
ing on the atomic transitions might provide an alterna¬ 
tive route to strong localization of light, as was recently 
suggested [55] . 

We acknowledge fruitful discussions with Nicola Pi- 
ovella, Tom Bienaime, Romain Bachelard, Guillaume 
Labeyrie, and Dominique Delande and the technical 
help from Louis Bellando, Ivor Kresic, Loi’e Lavenu, 
and Antoine Dussaux, and we thank Alain Aspect 
for his constructive comments on the manuscript. We 
also acknowledge financial support from the French 
Agence National pour la Recherche (project LOVE, 
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Fig. 3. Dependence of the subradiance time with the on-resonance optical depth bg and the laser detuning S. (a) 
Subradiant time constant Tsub (in units of Tat) as a function of bo = 3N/(koR)^ for different detunings of the probe 
laser. Almost all of the points collapse on a single curve, showing the linear scaling with bg. The dashed line is a linear 
fit Tsub/Tat = 1 + o:bo with the slope a as a free parameter. Excluding the data with bg > 120 we obtain a = 0.8. 
(b) From the same data, Tgub is plotted as a function of the detuning for three different values of bo, illustrating that 
Tsub is independent of the detuning, (c) Similarly, the subradiant relative amplitude Agub is plotted as a function of 
the detuning for the same three different values of bo- Subradiance is more important near resonance and decreases 
towards a plateau at large detuning. 
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Supplemental Material 


1. Experimental details 

A. Calibration of ho. 

As we cannot measure bo simultaneously to the data 
acquisition, the calibration is done in a separate time 
sequence, using absorption imaging, requiring a longer 
cycle without recapture of the atoms. We therefore in¬ 
crease the loading time of the MOT such that the fluo¬ 
rescence of the probe beam measured at large detuning 
is the same as in the measurement cycle, all other pa¬ 
rameters being unchanged. The absorption images then 
provide the transverse profile of the probe transmission 
T, which is related to the optical depth by T{x,y) = 
exp[—b'{x,y)], where 


b\S) = C 


bp 

1 -h 4(52 ’ 


( 1 ) 


with 6 the probe detuning and C = 7/15 is the average 
Clebsch-Gordan coefficient of the transition for a statisti¬ 
cal mixture of the Zeeman substates. We use a Gaussian 
fit for b'{x,y) and extract the r.m.s. size R directly from 
the fit, as well as bo from the amplitude of the Gaus¬ 
sian fit and the atom number N from its integral. This 
is done for different values for the time of flight of ballis¬ 
tic expansion, from which we extract the temperature, 
allowing us to know the size of the atomic cloud for the 
different pulses as well as the initial atom number. 

During the pulse series, we need to take into account 
possible optical pumping by the probe beam into the 
other hyperfine ground state F = 1. This effect is almost 
negligible near resonance but becomes significant for the 
largest negative detunings used in this work. We there¬ 
fore have precisely measured the pumping rate, from 
the fluorescence decay, at very large negative detuning 
((5 = —9), and used this rate to calibrate the probe in¬ 
tensity, which is measured simultaneously to the data 
on a separate photodetector. Then we use this intensity 
to compute the atom losses induced by optical pumping 
and correct the atom number for each pulse. 

Note that the uncertainty on the precise value C cor¬ 
responding to the experimental conditions introduces a 
systematic uncertainty on N and bo, which is thus also 
affecting the precise value of the slope extracted from 
the data of Fig. 3. Statistical uncertainties, estimated by 
the shot-to-shot fluctuations of bo, are of the order of 8% 
(standard deviation). 


B. Analysis of the fluorescence decay. 

The raw data are histograms of number of photons de¬ 
tected with a time bin of 1.6 ns. Our detection in the 
photon counting regime yields about 10^ counts per sec¬ 
ond (cps) during the pulses, while between pulses, after 
the end of the subradiant decay, the ‘dark’ level is about 


4000 cps, mainly due to stray light. For each pulse we 
subtract this dark level and normalize the signal to the 
steady state amplitude. Then we further bin the data 
to improve the signal to noise ratio. We use a variable 
binning, i.e., a larger bin size for later time, when the 
signal varies slowly, than for short times. We also imple¬ 
mented a systematic procedure to choose the range of 
data to be htted by selecting data one decade above the 
noise floor. This allows us to fit the last measurable time 
constant. We found this procedure to give the most reli¬ 
able and relevant fit. Finally, we compute the statistical 
coefficient of determination of the fit, which quantifies 
its relevance, and we keep only the points for which this 
coefficient is above 0.97 in Fig. 3 of the main paper. 

The data gathered in this figure have been taken in ex¬ 
actly the same conditions (initial atom number, tempera¬ 
ture, saturation parameter, experimental cycle), where 
the only changing parameters are the probe detuning 
and its intensity to keep the saturation parameter con¬ 
stant. We have also been able to observe subradiant de¬ 
cay for larger detunings, up to i5 = —11, but we have not 
included these data because we had to adapt other pa¬ 
rameters. Indeed, because of the limited available power 
for the probe beam, the saturation parameter and the 
subsequent amount of fluorescence were lower, which we 
partially compensated by increasing the atom number 
via the loading time, which changed also the tempera¬ 
ture and the cycle duration. 


2. Predictions of the coupled-dipole model 

We recall here the main ingredients of the coupled-dipole 
model, which has been widely used in the last years in the 
context of single-photon superradiance [iHin]- We con¬ 
sider N two-level atoms (positions r^, transition wave¬ 
length A = 2'K/ko, excited state lifetime F“^) driven 
by an incident laser (Rabi frequency fl, detuning A, 
plane wave of wavevector ko). Restricting the Hilbert 
space to the subspace spanned by the ground state of 
the atoms |G) = \g ■ ■ ■ g) and the singly-excited states 
|i) = \g ■ ■ ■ Ci ■ ■ • g) and tracing over the photon degrees 
of freedom, one obtains an effective Hamiltonian describ¬ 
ing the time evolution of the atomic wave function |^(t)), 


N 

m)) = oc{t)\G)+Y,m\^- ( 2 ) 

i=l 


Using standard approximations, the effective Hamilto¬ 
nian can be written as 


i7eff = — -f s\] 




( 3 ) 


i 
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where S'j- and SI are the usual pseudospin operators for 
the kets gi and Cj, respectively. The first term of 
describes the coupling to the laser field, the second ac¬ 
counts for the finite lifetime of the excited states and the 
third one describes the dipole-dipole interactions, with 


kra 


ri-i = \ri - Vi 


( 4 ) 


Here, we have considered a scalar model for light, which 
neglects polarization effects and near-field terms in the 
dipole-dipole interaction. It is known to be a good ap¬ 
proximation for dilute clouds [nun], i.e., when the typi¬ 
cal distance between atoms is much larger than the wave¬ 
length, which is the case in the experiment. 

Considering the low intensity limit, when atoms are 
mainly in their ground states, i.e., a ~ 1, the problem 
amounts to determine the amplitudes Pi, which are then 
given by the linear system of coupled equations 


'i= [iA- 


Pi - 


ir 


-J2v,,p,. (5) 

i¥=3 


These equations are the same as those describing N clas¬ 
sical dipoles driven by an oscillating electric field 0 , jus¬ 
tifying the term “coupled-dipole model”. The first term 
corresponds to the natural evolution of the dipoles (os¬ 
cillation and damping), the second one to the driving 
by the external laser, the last term corresponds to the 
dipole-dipole interaction and is responsible for all collec¬ 
tive effects, including dephasing and attenuation of the 
driving laser beam, as well as more subtle multiple scat¬ 
tering and cooperative effects (super- and sub-radiance). 
Note that even if the detuning does not appear explicitly 
in the dipole-dipole interaction term, it still strongly in¬ 
fluences the collective behaviour of the system through 
the population of the eigenmodes that contribute to the 
system response to the driving field. At large detuning 
for instance, multiple scattering vanishes. 

From the computed values of Pi, we can derive the 
intensity of the light radiated by the cloud as a function 
of time and of the angle [6] . The time dependence of the 
total radiated power P after switching off the laser is 
proportional to the derivative of the total excited state 
population, 

, N 

P ■ ( 6 ) 

We have used this model to study superradiance and 
subradiance decay in more detail than in the previous 
work [7]. The complete study will be published else¬ 
where m- Here, we only show the most important re¬ 
sult for the subradiance experiment, which is the lin¬ 
ear scaling of the decay time with the parameter 6o = 
3N/(koR)^, which also corresponds to the on-resonance 
optical depth through the center of the cloud. Here N is 
the atom number and R the r.m.s. radius of the atomic 
Gaussian density distribution, which are the two param¬ 
eters of the simulation. However, in order to determine 



Fig. 4. Scaling of the the subradiant decay time con¬ 
stant as a function of bg in the coupled dipole model. 
We used different densities po-\^ = {0.5,1, 3, 5, 7,9}. The 
time constant is extracted from the exponential fit of the 
end of the decay of the total emitted power computed 
from Eqs. (isrei) . The fit window is chosen for P (normal¬ 
ized to 1 at t = 0) to be between 10“'^ and 10“®. 


if bg is the scaling parameters, and not, e.g., the atomic 
density, we need to change bg while keeping the density 
p oc N/R^ constant, by changing simultaneously N and 
R. Since the atom number N is limited to a few thou¬ 
sands in our numerical simulations, quite large densities 
are required to reach large bg. It is then important to use 
an exclusion volume when drawing the random positions 
of the atoms in order to exclude close pairs of atoms [12] . 
We set our exclusion volume as kgVij > 3. 

In Fig. [4] we show the result of such simulations, for 
which we varied bg at constant density, for different 
values of pg\^ = N/{kgR)^, where pg is the peak 

density (at the center of the Gaussian cloud). We clearly 
see that the different curves collapse on the same line, 
showing that the on-resonance optical depth is indeed 
the parameter that controls the subradiant decay. The 
slope of the line is 0.36. Note that we do not a priori ex¬ 
pect this number to be in quantitative agreement with 
the experimental results, as polarization effects and the 
complex Zeeman structure of the rubidium atoms used 
in our experiment are neglected in the model. 

3. Test of different scaling parameters 

During the series of pulses, two effects contribute to a 
change of bg. The most important one is the expansion 
of the cloud due to its thermal velocity distribution. Its 
r.m.s. size typically varies from 0.5 to 1.2 mm between 
the first and the last pulse. This variation is independent 
of the probe laser detuning. Another contribution to the 
change of bg between subsequent pulses is the loss of 
atoms due to optical pumping into the F = \ ground 
state. Since we kept the same saturation parameter for 
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bo = 3N/{koRf N/{koR)^ 


Fig. 5. Test of different scaling parameters. The subradiant decay time constant is plotted as a function of different 
combinations of N and R. The collapsing of all points on a single curve is best with bo as the scaling parameter. 


all detuning, the number of scattered photons is always 
the same, about 25 per atom and per pulse. However, the 
probability that a scattering event transfers the atom 
to the F = 1 state depends strongly on the detuning 
from the F = 2 ^ F' = 2 transition, and thus it is 
much larger for large red detuning S (defined from the 
F = 2 F' = 3 transition). For example, at <5 = —6, the 
number of atoms in the F = 2 ground state has decreased 
by 70% between the first and the last pulse, whereas on 
resonance it is almost constant for all pulses. Thus, for 
the different detunings the relative contributions of the 
variation of iV and R are different. 

We have exploited this fact to test which combina¬ 
tion of N and R allows the different curves acquired 
with different detunings to collapse on a single curve. 
This is shown in Fig. [SJ where it clearly appears that 
the combination N/{koR)^ provides the best scaling pa¬ 
rameter. This experimentally confirms that bo is in¬ 
deed the parameter governing subradiance in dilute and 
extended samples, as predicted by the coupled-dipole 
model (Fig. |4]). 

4. Subradiance as a function of the probe inten¬ 
sity 

In order to validate the weak excitation assumption, we 
have also studied the subradiance decay as a function 
of the probe intensity. In the linear optics regime, the 
probe intensity should not impact on the measured de¬ 
cay rates or normalized amplitudes. We have thus fixed 
the detuning 6 and the optical depth and varied the in¬ 


tensity I. On Fig.|n]we plot the results for the subradiant 
time constant and relative amplitude as a function of the 
saturation parameter 


s(5) 


///sat 

l + 4S^ 


(7) 


where /gat is the saturation intensity of the transition. 

No significant dependence on the saturation parameter 
in the explored range (10“^ < s{6) < 10“^) is observed, 
confirming that the other experimental data, for which 
s(i5) ~ 0.05, have been taken in the linear regime. 


5. Possible role of radiation trapping of near¬ 
resonant light 

The very slow decays observed at large detuning cannot 
be simply interpreted as radiation trapping, as in the 
experiments by Labeyrie et al. [HI ITS] . since radiation 
trapping does not depend on bo only but depends on 
the optical depth ‘seen’ by the laser beam, i.e., b'{5) oc 
6o/(l + 4(5^), which is very small at large detuning. 

However, since the amplitude of the subradiance de¬ 
cay is very small, we have to insure that it cannot be 
explained by radiation trapping of a small amount of 
light that would always be on resonance. 

We have investigated two possible nonnegligible 
sources of near-resonant light: inelastic scattering due 
to the atomic saturation and the wings of the laser spec¬ 
trum. 









Fig. 6. Effect of the probe intensity. Subradiance time constant (a) and relative amplitude (b) as a function of the 
saturation parameter. The data have been taken with bo = 110 ± 8% and S = —5. 


A. Inelastic scattering 

The scattering rate Fsc of an atom driven by a laser is the 
sum of two contributions, elastic and inelastic scattering, 
given by m 

^ ^ f s 

bsc = rei+ri„ei = 2 (^(1 + 5)2 + )^2 TT7 ’ 

( 8 ) 

where s is the saturation parameter dehned above 
[Eq. d?!)]. Usually, inelastic scattering is neglected when 
s 1, which is the case in our experiment. However 
since subradiance is a small signal, special care is re¬ 
quired. Eor s ^ 1, the relative weight of inelastic scatte¬ 
ring compared to the total scattering is equal to s. The 
spectrum of this inelastic scattering is the well-known 
Mollow triplet which is, in fact, a triplet only at 
large s. At large detuning and moderate intensity, it is 
a doublet constituted of two lines of equal weights, one 
of which is precisely on the atomic resonance [18] . Thus, 
we expect that after the first scattering event, a propor¬ 
tion s/2 of the light is shifted near resonance. Multiple 
scattering could then trap this light in the sample. 

To estimate the weight of this inelastic component in 
the slow decay, we multiply this s/2 fraction by the rela¬ 
tive amplitude of the slow decay for resonant photon, 
which is what we extract from our experiment when the 
laser is tuned to resonance (Agub — 2 x 10“^ for bo = 110, 
see Eig. 3c of the main paper). In this scenario, we would 
thus expect Agub — s x 10“^ <C 10“^. On the contrary, 
we see on Eig. [SK also with bo = 110) that Agub is signifi¬ 
cantly larger and in any case not proportional to s. These 
data allow us to conclude that spurious resonant radia¬ 
tion trapping of one Mollow sideband is much smaller 
than the subradiance we observe. 

Note that in this experiment, a laser beam larger than 
the atomic cloud has been used, leading to a significant 
proportion of single and low-order scattering events on 
the edges of the cloud, even on resonance, contrary to 
the previous experiments on radiation trapping [CTIT^ . 


where a small beam at the center of the atomic cloud 
was used. That is why the relative amplitude of the slow 
decay is low even on resonance. 

B. Spectrum of the probe laser 

In our experiment, we used a distributed-feedback 
(DEB) laser diode as probe. This laser is expected to 
have a spectral linewidth on the order of ^ 3 MHz [19] 
corresponding to 0.5r, which should not have any in¬ 
fluence at large detuning. However, if the spectrum has 
slow-decaying wings (e.g., Lorentzian), the amount of 
light a few T’s from the central frequency might be not 
negligible in the experiment. 



Eig. 7. Spectrum of the probe laser. The power spectrum 
density (PSD) of a beat-note signal with a reference laser 
has been averaged 100 times (scanning time ~ 1 s) and 
the central frequency 1 GHz) has been shifted to zero 
(black dots). The result can be fitted by a Gaussian in 
the central part (red) and by a Lorentzian in the wings 
(blue). 
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Fig. 8 . Main results of a random walk model, (a) Temporal decay computed with t»'(0) = 53 and 5 = — 6 . The relative 
amplitude of the slow part is 10“"^. (b) Fitted time constant of the slow tail as a function of the on-resonance optical 
thickness. The increase is clearly faster than linear. 


To characterize the spectrum of our probe laser, we 
performed a beat-note experiment with a commercial 
extended-cavity laser diode (Toptica DL pro), with a 
much narrower spectral width. It is thus a good (but 
slightly conservative) approximation to consider the 
power spectrum density of the beat-note signal, recorded 
with an electronic spectrum analyzer, as the optical spec¬ 
trum of our probe laser. Such a measurement is shown 
in Fig. [T] The spectrum is composed of a central part, 
which can be well fitted by a Gaussian, superimposed 
on a small component with large wings, which is well 
fitted by a Lorentzian. The r.m.s. width of the Gaussian 
is 2.6 MHz (corresponding to a full width at half max¬ 
imum of IF, larger than expected) and the Lorentzian 
wings have a relative amplitude of 5.8 x 10“^ and a full 
width at half maximum 7 « 13.2 MHz corresponding to 

2 . 2 r. 

From the measured spectrum, we can numerically 
compute the amount of light near resonance (in a width 
F) given the detuning of the central frequency. For ex¬ 
ample, for 6 = —6r, corresponding to the largest detun¬ 
ing used in the measurements, this relative amount is 
ft! 1.5 X 10“"^. We recall that even when the laser is set 
on resonance, the relative amplitude of the slow decay is 
only Hsub ft 2 X 10“^ so that the two numbers should be 
multiplied to predict the relative amplitude of the slow 
decay induced by the small part of resonant light. This 
is thus much too small to be visible in the experiment. 

C. Random walk model with frequency redistribution 

One further spurious effect is frequency redistribution 
induced by the Doppler effect during the multiple scatte¬ 
ring process, which is known to play an important role in 
radiation trapping experiments in cold atoms [1^151120] . 
First, for light initially on-resonance, the light spectrum 
gets broader during the multiple scattering and the ef¬ 


fective optical thickness decreases, breaking the r oc 5^ 
scaling for high enough temperature. Second, light ini¬ 
tially slightly detuned has a tendency to drift towards 
resonance (in addition to the broadening). It means that 
the previous evaluation considering a window of width 
F to dehne the resonant light may be not appropriate 
because frequency redistribution induces a broader ‘cap¬ 
ture range’ of frequency. 

To evaluate these effects, we performed numerical sim¬ 
ulations considering the random walk of photons in a ho¬ 
mogeneous spherical cloud of atoms, taking into account 
the Doppler effect and its associated frequency redistri¬ 
bution. We send photons one by one, we randomly draw 
their initial transverse positions and detuning taking into 
account the spatial shape and measured frequency spec¬ 
trum of the probe laser, and we compute how many scat¬ 
tering events the photons undergo before escaping, which 
is converted to a time, using the fact that each scatte¬ 
ring event takes in average Tat [II]- We obtain thus a 
distribution of escape time, which is then convoluted by 
the duration of the pulse in order to simulate the fluo¬ 
rescence decay at the switch-off. 

In Fig.lHKa), we show the fluorescence decay for a laser 
detuning centered at 5 = — 6 , a cloud with a temperature 
T = 50 nK and an optical thickness b'{0) = 53, corre¬ 
sponding to bo ft 110 [Eq. (1)], for comparison with the 
data of Figs. 2-3 of the main paper. We see indeed that 
the small amount of light that is on resonance induces a 
slow decay, whose relative amplitude is ~ 10“^. As an¬ 
ticipated, this is larger than expected from the previous 
simple evaluation, but it is still more than one order of 
magnitude smaller than the subradiance decay observed 
in the experiment. Moreover, the use of a homogeneous 
sphere instead of a Gaussian density distribution overes¬ 
timates substantially this amplitude because for a given 
peak-optical thickness (as measured by the transmission 
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of a small beam at the center of the cloud), the average 
optical thickness ‘seen’ by a large beam is larger with 
a homogeneous sphere than with a Gaussian cloud. We 
have also checked that if we consider an incident short 
pulse of duration similar to the switch-off of the probe 
beam, and its associated spectral broadening, the result¬ 
ing slow decay also has a relative amplitude below 10 “^. 

Finally, we show in Fig. [8l)b) how the time constant 
associated with this small, slow decay, which can be fit¬ 
ted by a decaying exponential, evolves with the optical 
depth of the cloud. We can see that the increase with the 
optical depth is much faster than linear in this regime 
of low temperature [2^. This provides another, qualita¬ 
tive difference between our subradiance data and what 
could be expected from spurious radiation trapping of 
near-resonant light. 
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